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ABSTRACT 

Context. The origin and composition of Ultra-High Energy Cosmic Ray Events (UHECRs) are under debate. Possible sources include Active 
Galactic Nuclei - selected by various criteria - and extragalactic magnetars. 

Aims. We aim to improve constraints on the source population(s) and compositions of UHECRs by accounting for UHECR deflections within 
existing Galactic magnetic field models (GMFs). 

Methods. We used Monte Carlo simulations for UHECRs detected by the Pierre Auger Observatory and AGASA, to determine the UHECR 
trajectories within the galaxy and their outside-the-Galaxy arrival directions. The simulations, which used UHECR compositions from protons 
to Iron and seven models of the ordered GMF, accounted for uncertainties in the GMF and a turbulent magnetic field component. The trajectories 
and outside-the-Galaxy arrival directions were compared with Galactic and extragalactic sources. 

Results. For a given proton or light UHECR, the multiple potential outside-the-Galaxy arrival directions within a given GMF model are not 
very different, allowing meaningful constraints on source populations. The correlation between a subset of UHECRs and nearby extended 
radiogalaxies (Nagar & Matulich 2008) remains valid, even strengthened, within several GMF models. Both the nearest radiogalaxy Cen A, 
and the nearest radio-extended BL Lac, CGCG 413-019, are potentially sources of multiple UHECRs. The correlation appears to be linked to 
the presence of the extended radio source rather than a tracer of an underlying matter distribution. Several UHECRs have trajectories which 
pass close to or through the Galactic plane, some passing close to Galactic magnetars and/or microquasars. For heavier UHECRs, the multiple 
potential outside-the-Galaxy arrival directions of any given UHECR are highly scattered but still allow meaningful constraints. It is possible, 
but unlikely, that all UHECRs originate in the nearby radiogalaxy Cen A. 

Conclusions. Nearby radiogalaxies remain a strong potential source of a significant subset of UHECRs. For light UHECRs about a third of 
UHECRs can be "matched" to nearby galaxies with extended radio jets. The remaining UHECRs could also be explained as originating in 
extended radiogalaxies if one has at least one of: a large UHECR mean free path, a high cluster and/or intergalactic magnetic field, a heavy 
composition for two-thirds of the detected UHECRs. If extended radiogalaxies are, or trace, UHECR sources, the most consistent models for 
the ordered GMF are the BS-S and BS-A models; the GMF models of Sun et al. 2008 are acceptable if a dipole component is added. 

Key words, interstellar-medium: Cosmic Rays - interstellar-medium: Magnetic Fields - Galaxies: Active - Galaxies: jets; 



1. Introduction 

Ultra-High Energy Cosmic Rays (UHECRs) are protons or 
fully ionized nuclei with energies greater than about 10 '"^ 
eV (10 EeV). On entering the earth's atmosphere they pro- 
duce a shower of secondary particles and excite atmospheric 
molecules. The detection of both effects from the ground al- 
lows a precise measurement of both the initial energy and 
the arrival direction of the UHECR above the earth's at- 

1999)1 and the Pierre 



mosphere e.g. AGASA (Takeda et al. 



Auger Observatory ( Pierre Auger Collaboration et al. 2004[ l 



The resultant positional accuracy of the earth-arrival direction. 
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now better than a degree, allows correlations with astronomi- 
cal sources in order to determine the source population(s) of 
UHECRs - a long held mystery. Exploiting the full potential 
of the earth-arrival directions nevertheless requires taking into 
account deflections suffered due to Galactic and extragalactic 
magnetic fields. 

Suspected sources of ultra high energy cosmic rays include 
powerful active galaxies, magnetars, core collapse Supernovae 
(SN II), and gamma rays bursts (for reviews and lectures 
see |HiIlas| ( |1999| ); [KacheMess] ( |2008| l and references therein). 
Specifically, an origin in nearby radiogalaxies and BL Lacs has 



long been predicted (Rachen & Biermann 1993 Romero et 



al. 1996 1. More exotic explanations include dark matter anhi- 
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lation and evaporation of micro black holes. The main energy 
loss for a UHECR propagating over cosmological distances is 
expected to be pion-production triggered by interaction with a 
CMB photon, the so called Greisen Zatsepin Kuzmin (GZK) 



of pulsars and extragalactic radio sources, and polarization of 



starUght or radio synchrotron emission (e.g. Heiles 



effect (Greisen 1966 Zatsepin & Kuzmin 1966| l. The pre- 
dicted energy loss (around 30%) and mean free path of the 
UHECR (~20-100 Mpc) depend on the composition and en- 
ergy of the UHECR and details of the collision process, e.g. the 
photo-proton cross section (see e.g. Stanev ( 2004[ l; Kachelriess 
( |2008l l; PA08). 

The Pierre Auger Observatory (PAO) has detected 81 
events between 2004 January 1 and 2007 August 3 1 with recon- 
structed energies above 40 EeV and zenith angles smaller than 
60°; of these, 27 have energies above 56 EeV and have been re- 



ported in Pierre Auger Collaboration et al. ( 2007 2008 here- 
after PA07, PA08). The origins of the latter 27 UHECRs, to- 
gether with 1 1 UHECRs detected by AGAS A at energies above 
56 EeV (Hayashida et al. 2000| l are the focus of this paper. 
The PAO has since detected additional UHECRs ( |The Pie^e| 
Auger Collaboration et aT] |2009a|b| hereafter PA09a, PA09b) 
but detailed positions and energies are not yet published. In 
this work we use "UHECR" to refer to UHECRs with energy 
above 56 EeV unless explicitly mentioned otherwise. The com- 
positions of the detected UHECRs are still under debate (e.g. 



PA09a; PA09b; |Matthews|[2007t [Hooper & Taylor|[2009| ) with 
evidence for a mix of protons and heavy nuclei at the high- 
est energies. There has been recent doubt about the measured 
energies of the UHECRs detected by PAO and AGASA. The 
PAO team have lowered their previous reported energies by 
a small percentage (PA09a, PA09b) and there have been sug- 
gestions that energies reported by AGASA should be lowered 
by up to 30%. To avoid confusion we report our results sepa- 
rately for PAO and AGASA UHECRs. Several other UHECR 
observatories, e.g. Haverah Park and Fly's Eye, have detected a 
significant number of UHECRs: given that the GMF deflection 
calculations and source matching in this work require accurate 
energies (<10%) and arrival directions (<2°) we do not include 
their results in our analysis. 

The arrival directions of the 27 UHECRs detected by PAO 
are not isotropic at a 99% significance level (PA07, PA08, 
PA09a, PA09b). Several correlations, or lack thereof, with ex- 
tragalactic sources have been suggested: AGNs from the cata- 
log of Veron-Cetty & Veron (PA07, PA08), hard X-ray selected 
AGNs ([G eorge et al. 2008), nearby extended radiogalaxies 
( Nagar & Matulich,, 2008 j ), nearby spiral galaxies (Ghisel lini et| 
al. 2008 1, local volume galaxies (Cuesta & Prada 2009) and 
large scale structure ( Ryu et al. 2009). 

Galactic and extragalactic magnetic fields deflect cosmic 
rays via the Lorentz force even at energies above 56 EeV. If the 
extragalactic field is unordered over large scales the expected 
deflections for these energies are only a few degrees (e.g. 
Harari et aL) [2002[ [Dolag et al] [2004| l. An ordered Galactic 
Magnetic Field (GMF) can, however, produce significant de- 



flections (e.g. Takami & Sato 2008 1 especially on trajectories 
which pass close to the Galactic plane and/or for heavy com- 
position UHECRs. 

The line of sight integrated GMF of the Galaxy has been 
measured in various directions via e.g. the rotation measures 



[Brown et al.] [2007i 'Noutsos et al] [20081 [Sun et al. 





1996 




2008 



Han 2009 Beck 2009). Additionally, local magnetic fields 



can be measured accurately via Zeeman splitting of emission 
lines. For a given line of sight, the rotation measure traces the 
integrated parallel component of the GMF while a UHECR 
is deflected by the perpendicular component of the GMF. A 
global model of the ordered GMF is thus required to trans- 
fer information from rotation measures into information useful 
for UHECR deflections. Several models have been presented 
for ordered magnetic fields in the Galaxy. These typically con- 
sider one or more of the following: a toroidal type field in and 
near the galaxy disk, a toroidal field in the halo, and a poloidal 
field which reproduces the vertical component seen in the Solar 
neighbourhood and the Galactic center. Models for the ordered 
component of the magnetic field are discussed in the next sec- 



tion, and recent comprehensive reviews can be found in Han 
( [2008| and [Beck | (pons). Studies of the global magnetic field 
of other galaxies has greatly faciUtated the refinement of mod- 
els of the Galactic GMF (see Beck) [20081 |2009l for a review), 
even though there is significant variation between galaxies in 
both the ordered and turbulent GMF, especially among those 
with high star formation rates. 



In this article, we expand the work of Nagar & Matulich 
(j2008| in two directions: first we use Monte Carlo simulations 



to derive UHECR trajectories within the Galaxy and corre- 
sponding outside-the-Galaxy arrival directions for all PAO and 
AGASA UHECRs. The simulations, mn for six UHECR com- 
positions - from protons to iron - and seven GMF models, in- 
clude uncertainties in the GMF and a turbulent magnetic field 
component. We then compare the UHECR trajectories within 
the galaxy and their outside-the-Galaxy arrival directions with 
various Galactic and extragalactic sources, and discuss the im- 
plications of the results. The individual steps followed in this 
work have been addressed previously by various authors ref- 
erenced here and in the next section. The new facet of this 
work is that we combine all of the above steps within the 
same Monte Carlo simulation, in order to model the trajecto- 
ries of the PAO and AGASA UHECRs with energies gi-eater 
than 56 EeV within a large number of ordered GMF models 
and UHECR compositions. Sec.|2]introduces the magnetic field 
models and describes the Monte Carlo simulations, Sec.[3]sum- 
marizes the sources of the data used, and Sec. [4] describes the 
principle results obtained. Finally, Sec. |5] contains a brief dis- 
cussion and the conclusions of our study. Distances to galaxies 
are calculated using a Hubble constant of 72 km s"' Mpc"', 
except for relatively nearby galaxies for which we use distances 
as referenced. 

2. Magnetic Fields and IVIonte Carlo Simulations 

We consider seven ordered magnetic field models for our 
Galaxy (for recent reviews of Galactic and extragalactic or- 
dered magnetic fields see [Beckl [20081 [Haiil [2008 : Beck| [2009 



Four of these models, the so called BS-S, BS-A, AS-S, and 
AS -A have been used previously by several authors. Here, the 
first two alphabets signify a (bi-)symmetric (BS) or asymmet- 
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ric (AS) configuration w.r.t. a transformation between and 
O + n, where is the cylindrical azimuthal angle in the Galactic 
plane. The last alphabet signifies whether the field reverses di- 
rection (A) or not (B) in passing from above the Galactic plane 
to below the Galactic plane. Given the evidence for a vertical 
(z) component of the GMF, both in the Galactic center and in 
the Solar vicinity, a dipole field model is usually added to the 
toroidal field above. Detailed descriptions and discussions of 
these GMFs can be found in 'Stanev'(f l^97]l;[Becp(|2 00B;'Han' 
(|2002|l;|Tinyakov & Tk achev (2002); iProuza & Smida (2003 ); 
Brown et al.| (|2007 ); KacheMeB e t aLl ( |2007| l; |Men et al.| f2008 ); 



Takami & SatO| ( i2008| ); |Han| ( |20(j9l l. For these four models, we 



use the equations and their related normalizations strictly fol- 



lowing Takami & Sato ( 2008 1, and always include the dipole 



field component as specified in Takami & Sato (2008 1. 

Sun et al.| (|2008| hereafter S08) have recently presented 
three new variations of the above GMF models - AS-S-i-RING, 
AS-Sh-ARM and a modified BS-S - tailored to best fit observa- 
tions of Galactic synchrotron emission and its polarization, and 
rotation measures of pulsars. For brevity we sometimes refer to 
these three models as ASS-hR, ASS-hA, and BSS(S) in the text 
and tables. [Sun et al.| ( [2008| ) find that the first two models pro- 
vide the best fit to current radio observations. These two mod- 
els use the basic configuration of the AS-S model, but introduce 
field reversals in a specified ring or arm, respectively. The mod- 
ified BS-S model is based on the traditional BS-S model but 
with different parameters and normalizations. The three new 
models presented by S08 also include a toroi dal halo field as 
initially presented in Prouza & Smida ( 2003 1, but do not in- 



clude a dipole magnetic field component. We use these three 
models and their related halo field as described in S08 with only 
one difference: in the AS-S+ARM model we directly use the 
log spiral arm morphology described in |Wainscoat et al.| ( |1992| 
without correcting for their shape near the solar neighbour- 



hood as in S08, who used the correction suggested in Taylor 



& Cordes ( 1993 i. We do not expect significant changes in the 



UHECR trajectories due to this small change in the morphol- 
ogy of the spiral arms. [Sun et al.|p008[ l assumed that there is 
no ordered dipole field component in the Galaxy GMF, argu- 
ing that the vertical (z) component of the local magnetic field 
is a turbulent component. Since the parametric form and nor- 
malizations of the S08 models were derived via fits to radio 
observations, it is not strictly correct to add a dipole field to the 
models. Nevertheless, for illustrative purposes (and out of cu- 
riosity) we also performed simulations of the S08 GMF models 
including the dipole field component as specified in [Takami & 



[Sato] poos] ). We refer to these results as the Sun et al. (2008) 
GMFs plus dipole. 

We simulated UHECR trajectories for six different UHECR 
compositions: protons, He^^, O^^, Al'^*^, Ca^°^, and Fe^''^. 
Once chosen, a composition was maintained fixed during the 
full trajectory within the galaxy: i.e. we do not consider photo- 
dissociation-produced composition changes to a UHECR dur- 
ing its Galactic trajectory. UHECR trajectories were simulated 
by firing an anti-particle with the corresponding UHECR mass, 
(anti)charge, and energy from the Earth in the earth-arrival di- 
rection of the UHECR and following its trajectory - including 
deflections in the magnetic field - in steps of 10 pc. The simu- 



lations were carried out until the anti-particle left a sphere with 
diameter 40 kpc centered on the Galactic center The Galactic 
coordinate pair (I, b) derived from the final position and veloc- 
ity of the anti-particle in the simulation is equivalent to the ar- 
rival direction of the UHECR before it entered the Galaxy, and 
we refer to this direction as the 'outside-the-Galaxy' arrival di- 
rection, to distinguish it from the 'earth-arrival' direction mea- 
sured by the observatory. 

Each UHECR's trajectory within each GMF model was 
simulated fifty times using a Monte Carlo approach: introduc- 
ing uncertainties in the ordered magnetic field and a random 
turbulent magnetic field component. In more detail, in each 
10 pc simulation step we (a) computed the 3-vector repre- 
senting the Cartesian coordinates of the ordered magnetic field 
from the relevant model. For the AS-S, AS-A, BS-S, and BS- 
A models this was the sum of the cylindrical (disk) field and 
the dipole field. For the AS-S+ARM, AS-S-i-RING, and BS- 
S model of S08 this was the sum of the disk and halo fields 
(with or without the dipole field as described above); (b) added 
a Gaussian-distributed random error in each component of the 
ordered magnetic field, i.e. a change in both magnitude and di- 
rection of the ordered field. This was done by adding to the 
3-vector of the previous step a 3-vector whose each component 
was a random number derived from a Gaussian distribution 
with mean = and cr - 0.5 times the value of the corresponding 
Cartesian component of the GMF; (c) added a turbulent field 
component. The turbulent fiel d component w as computed fol- 
lowing the recipe described in|Prouza & Smida (2003), though 
with slightly different normalization parameters: the maximum 
turbulent magnetic field strength in each Cartesian component 
was taken to be ±10/zG, and the probability of a turbulent field 
being present was taken to be 1 % in the halo of the galaxy, 20% 
in the disk of the galaxy (Galactic radius less then 20 kpc and 
Galactic height less than 1.5 kpc) and 80% in the spiral arms 
of the Galaxy - where we used the spiral arm definition from 



Wainscoat et al. (1992 1, spiral arm thickness 0.75 kpc in the 
disk (i.e. coordinate r), and half-height 0.5 kpc (coordinate z). 
We therefore generated a random number to determine whether 
a turbulent field was required to be added following the prob- 
abilities above. If so, we generated a 3-vector in which each 
component was a random number uniformly distributed be- 
tween ilOyuG. This was added to the magnetic field 3-vector in 
the previous step in order to derive the total magnetic field used 
in the deflection calculation of this step. This turbulent field 
was then kept constant for the next 50 pc (i.e. an additional five 
simulation steps), the expected scale length of tu rbulent field in 



the Galaxy ( Prouza & Srmda 



2003 



Beck 



2008 1. Note that the 



turbulent field used here is relatively large in comparison with 
current measurements of the turbulent field in the Galaxy (S08) 
and other galaxies ( [Beck [2009[ ). 



In summary, within each GMF model we derived 50 
outside-the-Galaxy arrival directions for each UHECR of a spe- 
cific composition, or a total of 300 outside-the-Galaxy arrival 
directions for each UHECR over all six compositions mod- 
elled. 
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3. Data 

Positions and energies of the 27 UHECRs detected by the PAO 
with energies > 56 EeV were taken from PA08; we have not 
made the small energy corrections to these UHECRs as later 
reported by the PAO (PAO09b). Positions and energies of the 
UHECRs detected by AGASA with energies > 56 EeV were 
taken from Hayashida et al. ('2000') which is an updation of 
the data listed in Takeda et al. (1999). For AGASA UHECRs, 



we lowered the energies reported in Takeda et al. (1999 1 by 
30%: this resulted in 11 UHECRs with energy above 56 EeV. 
In this work we use 'AGASA UHECRs' to refer to only these 
1 1 UHECRs unless explicitly stated otherwise. 

Our catalogs of astronomical sources were drawn from var- 
ious publications and web-based source lists, and chosen to be 
representative of previously suspected or posited Galactic and 
extragalactic sources of UHECRs. In all cases, we attempted 
to select the most comprehensive catalog available. Source 
lists for galaxies, galaxy clusters, GRBs, Galactic supernova 
remnants, and extragalactic radio supernova are described in 
Nagar & Matulich| ( ,2008) . As in |Nagar & Matulich] ( [2008] l, in 



the case of radiogalaxies and radio sources we used several 
surveys and catalogs, including NVSS ( [Condon et al.| |1998| l, 
SUMMS ( |Bock etal.|[T999i ), and NED. The list of extragalac- 
tic jets is based on that compiled by |Liu & Zhang (2002); to 
this list we added data on the galaxies as given in NED, and 
re-measured the total flux and total extent of the radio emission 
from NVSS and SUMMS maps or from the Ust of DRAGNs 
compiled by P. Leahy The DRAGNs list is a subset of the 



3CRR sources (Lamg et al. 1983) which have one or all of 



radio-emitting jets, lobes, and hotspots. The total extent of the 
radio emission therefore includes both radio jets and any ra- 
dio lobes. This extent is typically referred to as the Largest 
Angular Size (LAS) or Largest Linear Size (LLS). In the case 
of double sided jet or lobe sources, we added the LLS of the 
two jets and lobes as scalars instead of vectors in order to dis- 



count the effects of jet and lobe bending. Nagar & Matulich 



( 2008| ) listed the 10 radiogalaxies in the field of view of PAO 
with D< 75 Mpc and LLS > 180 kpc. To this we added all other 
(northern) radiogalaxies satisfying the same criteria: NGC 315 
(D= 69.5 Mpc; LLS= 1111 kpc), NGC 383 (D= 71.7 Mpc; 
LLS= 907 kpc), NGC 1275 (D= 74.1 Mpc; LLS= 468 kpc), 
NGC 5127 (D= 68.3 Mpc; LLS= 230 kpc), and CGCG 514- 
050 (D= 72 Mpc; LLS= 506 kpc). Together, these 15 galaxies, 
form our sample of all nearby (D< 75 Mpc) extended (LLS 
> 180 kpc) radiogalaxies. In the case of nearby galaxies, we 
also used the HIPASS catalog (IMeyer et al.', '2004[ [Wong et 



al. 2006 ) as used by Ghisellini et al. (2008) in their study of 



the correlation between UHECR arrival directions and nearby 
galaxies. Additionally, we used the revised Third Reference 
Catalog of Bright Galaxies (Corwin et al. 1994| ) to select a 
sample of nearby elliptical galaxies. 

In this work we have considered several new source 
catalogs: Galactic Soft Gamma-Ray Repeaters (SGRs) and 
Galactic Anamolous X-ray Pulsars (AXPs) were taken from 
[Woods & Thompson] ( |2006[ ) and updated with the online list 

' http://www.jb.man.ac.uk/atlas/dragns.html 



maintained by the McGill pulsar grouf|^ As of June 2009, this 
list includes six SGRs (two of which are candidates) and 10 
AXPs (one of which is a candidate). These SGRs and AXPs 
are believed to be high magnetic field radio pulsars or "mag- 
netars". Of the total of 16 SGRs and AXPs, 14 are close to the 
plane of the Galaxy, and two are in the SMC and LMC. 

Confirmed Galactic microquasars were taken from the list 



in Paredes ( 2005 ), and candidate micro-quasars from Combi 



|et al.| ( |2008 1. Various gamma-ray catalogs have been used: 
the BeppoSAX catalog of GRB and X-ray afterglows ([de| 



Pasquale et al. 2006 1, the BeppoSAX complete catalogue of 
GRBs (|Vetere et al. 2007| , Gamma-ray Blazars in northern sky 
( Sowards-Emmerd et al.||2003 ), Blazar counterparts for 3EG 
sources (|Sowards-Emmerd et al. 2004[), G amma-ray Blazar 



candidates (Sowards-Emmerd et a 



2005[ ), Sources detected 



by I SGRI (Bodaghee et al. 2007] ), and the Third EGRET cata 
log ([Hartman et al. |1999| ). We also used the HESS ( [Hofmann 



|2005| catalog of 54 Gamma-ray sources between 100 GeV and 
100 TeV as obtained from the HESS online catalog in June 
200SE1 



4. Results 

We performed several Monte Carlo test runs varying the mag- 
nitudes of the error of the ordered field and the maximum tur- 
bulent field strength. Here we present and discuss the results of 
simulations which used a Gaussian-distributed (mean - and 
cr= 50% of the value of the Cartesian GMF component) error 
in each Cartesian component of the ordered GMF, and a tur- 
bulent magnetic field uniformly distributed between ilOyuG in 
each Cartesian component of the magnetic field, with details as 
described in Sec. |2] We note that for the simulated GMF error, 
using 25%-100% instead of 50% as described above does not 
significantly change the trajectories in the case of anti-protons. 
The relative deflections suffered by the UHECR from the or- 
dered field and turbulent field components can thus be roughly 
judged by comparing the average and rm.s. of the deflection 
angles corresponding to the fifty outside-the-Galaxy arrival di- 
rections for each UHECR within each GMF model. 

We used the following distance limits between a single 
outside-the-Galaxy arrival direction and an extragalactic astro- 
nomical object in order to be considered a "match": 2.5° for H, 
3° for He, 6° for O, 9° for Al, 12° for Ca, and 15° for Fe. These 
semi-arbitrary values were chosen considering a constant addi- 
tive error of 2° - to account for e.g. the error in determining the 
arrival direction of the UHECR - plus a composition-weighted 
error of 0.5 times the charge of the UHECR which roughly ac- 
counts for errors in the deflection calculations and deflections 
by intergalactic magnetic fields but avoids too large "match 
radii" for heavy nuclei. To match an astronomical source with 
a UHECR we require that at least 4 of the 50 outside-the- 
Galaxy arrival directions are 'matched' by the above criterion; 
this >8% match probability was chosen as it implies roughly 
that a match cannot be ruled out at better than a 2cr level within 
the simulation parameters. 



^ http://www.physics.mcgill.ca/ pulsar/magnetar/main.html 
' http://www.mpi-hd.mpg.de/hfm/HESS/pages/home/sources/ 
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Graphical representations of our simulation results are 
shown in Figures [T]to[5] Coordinates of the arrival directions of 
the UHECRs detected by PAO (blue for UHECRs with energy 
> 75 EeV and green for UHECRs with energy between 56 and 
75 EeV) and AGASA (red) are shown as open circles of (ar- 
bitrary) radius 3.5°. The outside-the-Galaxy arrival directions 
of the UHECRs are plotted with small colored dots (50 per 
UHECR, each representing the result of one Monte Carlo sim- 
ulation). Figures [T] and [3] appear in the printed version of this 
manuscript while Figures |2] |4] and |5] are available in the elec- 
tronic version only. The outside-the-Galaxy arrival directions 
were plotted in the order red, green, blue, and yellow. When 
one of these colors is not seen for a UHECR, it lies below one 
of the following colors in the sequence above. For example, the 
pairs BS-S (red), BS-A (blue) and AS-S (green), AS-A (yel- 
low) have the same deflection for positive values of Galactic 
latitude (b). To avoid overcrowding, results for the three mod- 
els of SOS are shown separately from the other four models. 

Figures[T]to|5]include the positions of all galaxies with ra- 
dio jets from the catalog of Liu & Zhang (2002 1 within 500 kpc 



and the 16 confirmed and candidate SGRs and AXPs in the 
Galaxy ( |Woods & Thompson, 2006) . As in [Nagar & Matulich] 
( 2008| l - but with different cutoff^s - we have divided the galax- 
ies with radio jets into three redshift bins: "nearby" (D < 
75 Mpc; red circular symbols), "intermediate" (75 Mpc < D < 
200 Mpc; blue circular symbols), and "distant" (200 Mpc < D 
< 500 Mpc; black circular symbols). In all redshift bins, we dis- 
tinguish between galaxies with radio structures more extended 
than 180 kpc ("extended"; solid symbols) and those with radio 
structures less extended than 180 kpc ("compact"; open sym- 
bols). Filled triangles are used for Galactic magnetars: red for 
SGRs and blue for AXPs. Of the 16 magnetars, two are in the 
LMC and SMC and the others are close to the plane of the 
Galaxy. 



4.1. UHECR Trajectory Simulations: Protons 

Columns 5 to 7 of Tables [T] and |2j present the results of the de- 
flection simulations for proton UHECRs within the BS-A GMF 
model. Tabular results for other models are not shown for space 
reasons and can be requested from the authors. The first four 
columns of the tables identify the UHECR by an index num- 
ber, its arrival direction in Galactic coordinates, and its energy. 
Then, for each UHECR composition simulated, we list the av- 
erage and r.m.s. of the deflection angles corresponding to the 
50 outside-the-Galaxy arrival directions of each UHECR, fol- 
lowed by the nearby extended radiogalaxy which best matches 
the fifty outside-the-Galaxy arrival directions (see the previous 
section for a definition of a match). A nearby extended radio- 
galaxy is only listed if more than four of the 50 outside-the- 
Galaxy arrival directions are matched. For proton UHECRs, the 
average deflection and its r.m.s. for each UHECR are both typ- 
ically l°-5° though values greater than 30° are seen for a few 
UHECRs whose trajectories pass close to the Galactic plane 
and/or Galactic center. 

We emphasize four important results of the proton-UHECR 
trajectory simulations: (a) even after considering errors in the 



model GMF and adding a relatively strong turbulent mag- 
netic field component, the outside-the-Galaxy arrival directions 
of the majority of proton-UHECRs are typically sufficiently 
concentrated to allow meaningful comparisons with Galactic 
and extragalactic sources; (b) for a given proton-UHECR, the 
outside-the-Galaxy arrival directions are typically significantly 
diff'erent for diff'erent GMF models. An assumption of the 
source population therefore allows constraints on competing 
models of the ordered GMF; (c) almost all UHECRs with 
Galactic latitude -25°< b < 0°have trajectories which pass 
close to or cross the Galactic plane within certain GMF mod- 
els, especially those which include an ordered dipole field; 
(d) as discussed in detail in Sec. |4.5| the match between 



proton-UHECRs and nearby extended radiogalaxies ( [Nagar &" 
Matulich 2008|) is maintained in most GMF models, and sig- 
nificantly strengthened in some GMF models. 



4.2. UHECR Trajectory Simulations: Heavy Nuclei 

With increasingly heavy UHECR compositions, the increased 
deflections suff'ered by the UHECRs tend to scatter their out- 
of-Galaxy arrival directions over larger regions of the sky. 
Figure[3] which compares the deflections suffered by UHECRs 
in four models for different UHECR compositions, empha- 
sizes the importance of UHECR composition for source iden- 
tification. The simulation in this figure does not include a 
Monte Carlo approach: i.e. we did not consider GMF uncer- 
tainties or turbulent magnetic fields. We remark on several in- 
teresting features in Figure [3] especially for compositions near 
Oxygen: (a) the three PAO UHECRs near I = 200°, b = -45° 
move close to the radiogalaxy 3C120 for a UHECR compo- 
sition of Oxygen; (b) the twin mystery PAO UHECRs near 
I = 335°, b - -18° move toward the Galactic plane in the 
case of light compositions, and then toward the Cen A concen- 
tration of radiogalaxies for an Oxygen composition, (c) several 
UHECRs move toward the extended radiogalaxies NGC 315 
and NGC 383 (near position / = 130°, b = -30°) for composi- 
tions near Oxygen. 

A summary of the results of our Monte Carlo simulations 
for heavy-nuclei UHECRs within the BS-A GMF model can be 
found in Tables [T] and |2] Graphical representations for Oxygen 
and h-on UHECRs within the BS-S, BS-A, AS-S, and AS-A 
models are presented in Figures |4] and |5] As expected, both the 
average and the r.m.s. of the deflection angles corresponding to 
the 50 outside-the-Galaxy arrival directions for each UHECR 
increase with increasing UHECR mass (Tables [T] and |2]i. For 
compositions up to around Oxygen, the 50 outside-the-Galaxy 
arrival directions corresponding to each UHECR are typically 
still identifiably together on the sky - allowing potential source 
searches (Figure|4]l. For compositions close to and higher than 
Ca, the fifty outside-the-Galaxy arrival directions correspond- 
ing to each UHECR are typically scattered over most of the sky. 
However, even in the case of Iron UHECRs (Figure |5]l a con- 
centration can still be seen in the region of the SuperGalactic 
plane near to and north of Cen A. 
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4.3. Source Populations: Centaurus A 

The concentration of UHECRs in the Cen A region has drawn 
a lot of attention since its first report by the PAO team. Our 
Monte Carlo simulations can be used to test whether Cen A is 
the dominant, or only, source of UHECRs. The results of these 
tests are shown in Table [3] for UHECRs detected by PAO and 
Table |4]for UHECRs detected by AGASA. In these tables, the 
first column for each composition reports the fraction of the 
total (1350 for PAO and 550 for AGASA) outside-the-Galaxy 
arrival directions which fall within three 20° circles centered on 
the Cen A nucleus and its northern and southern radio lobes. 
The second column for each composition shows the number 
of UHECRs (of the total of 27 for PAO and 11 for AGASA), 
for which at least one of the fifty outside-the-Galaxy arrival 
directions falls in the 20° circles described above. 

Within our simulation parameters, it is possible but unlikely 
that Cen A is the source of all PAO-detected UHECRs: for pro- 
ton UHECRs one third or less of the total outside-the-Galaxy 
arrival directions are within the 20° circles around Cen A in 
any GMF model, and almost two thirds of PAO UHECRs do 
not have even one outside-the-Galaxy arrival direction within 
the above area. For heavier compositions a smaller fraction of 
total outside-the-Galaxy arrival directions are within the 20° 
circles around Cen A but a larger fraction of PAO UHECRs, al- 
most 100% for Ca and Fe, have at least one outside-the-Galaxy 
arrival direction which falls within the Cen A circles. The lat- 
ter result is due to the large scatter of the outside-the-Galaxy ar- 
rival directions of any given heavy composition UHECR. In the 
case of AGASA UHECRs, the fraction of outside-the-Galaxy 
arrival directions which fall within the Cen A circles is very 
small for all compositions, though for heavy UHECR composi- 
tions, a large fraction or all UHECRs have at least one outside- 
the-Galaxy arrival direction falling within the Cen A circles. 

If we consider all 300 outside-the-Galaxy arrival directions 
corresponding to each UHECR (6 compositions, 50 outside- 
the-Galaxy arrival directions per composition) then within 
the models AS-S, AS-A, BS-S and BS-A, all PAO-detected 
UHECRs have at least one outside-the-Galaxy arrival direction 
which falls within the 20° circles around Cen A. In the Sun et 
al. 2008 models a few (two to six) PAO UHECRs do not sat- 
isfy the above criterion. For the AGASA-detected UHECRs, 
for the AS-S and AS-A models all UHECRs have at least one 
outside-the-Galaxy arrival direction which falls within the 20° 
circles around Cen A for some one composition, while in the 
other models between two and eight UHECRs do not satisfy 
this criterion. 



4.4. Source Populations: The Radio Galaxy 
Concentration around Cen A 



As Nagar & Matulich ( 2008[ l have pointed out, the area around 
Cen A hosts the highest density of nearby extended radiogalax- 
ies. We repeated the analysis of the previous section, but using a 
box around all nearby extended radiogalaxies in this area. For 
simplicity we used the following limits for the box: Galactic 
longitude between 240° and 360° and Galactic latitude be- 
tween -10° and 70° (the "RG box"). The results are hsted in 



sub-columns 3 and 4 for each composition in Table [3]for PAO 
UHECRs and Table |4]for UHECRs detected by AGASA. 

As expected the match statistics are better than in the case 
of only the 20° circles around Cen A and its radio lobes. The 
"RG box" contains between a quarter and a half of all outside- 
the-Galaxy arrival directions for any one given composition. 
For light UHECRs at least half of all PAO-detected UHECRs 
have at least one outside-the-Galaxy arrival direction which 
falls in the RG box. For Ca and Fe compositions all UHECRs 
satisfy this criteria. In the case of AGASA UHECRs heavy 
compositions are required to provide a good match between 
the UHECRs and the RG box. 

If we consider all 300 outside-the-Galaxy arrival directions 
corresponding to each UHECR then within the models AS-S, 
AS-A, BS-S and BS-A, all PAO and AGASA UHECRs have 
at least one outside-the-Galaxy arrival direction which falls 
within the RG box. In the Sun et al. models a few (two to four) 
AGASA-detected UHECRs do not satisfy the above criterion. 

4.5. Source Populations: Nearby Extended Radio 
Galaxies 



|Nagar & Matulich] ( |2008| l have previously argued that the ar- 
rival directions of a subset of PAO-detected UHECRs are cor- 
related with the directions of nearby extended radiogalaxies. A 
distance of less than 3.5° was found between the radio struc- 
tures of six extended radiogalaxies (of the 10 "visible" to the 
PAO) and 8 UHECRs detected by PAO. Considering the full 
sky, we have a total of 15 nearby (D< 75 Mpc) galaxies with 
extended (> 180 kpc) radio structures, with the new galax- 
ies listed in Sec. [3] to be compared to 38 UHECRs. Adding 
the AGASA-detected events, adds two new matches to earth- 
arrival directions of UHECRs: an AGASA event with original 
reported energy E= 120 EeV matched to NGC 7626 (which is 
also matched to a PAO-detected UHECR) and a match between 
an AGASA event with original reported energy E= 68 EeV 
with CGCG 514-050 (this event does not appear in the fig- 
ures as its 30% reduced energy is below our UHECR cutoff 
of 56 EeV). The correlation between nearby extended radio- 
galaxies and UHECRs detected by PAO and AGASA therefore 
remains highly statistically significant even before consider- 
ing deflections by the GMF. We remark that the Supergalactic 
plane in the Cen A region passes close to the Galactic longi- 
tudes where the deflection produced by an ordered GMF with 



azimuthal symmetry is minimal (e.g. Takami & Sato 2008 1. It 
may be this fortuitous coincidence which allowed us to find 
several UHECR events "matched" to nearby extended radio 
galaxies in [Nagar & Matulich| ( |2008| l. 

Including deflections of UHECRs by the GMF main- 
tains, and in some models increases, the correlation between 
UHECRs and nearby extended radiogalaxies. Table |5] in its 
first four columns, lists the statistics of matches between 
UHECRs and our sample of 15 nearby extended radiogalax- 
ies for all compositions and models simulated. The BS-S and 
BS-A models clearly result in the best match between nearby 
extended radiogalaxies and UHECRs. For proton UHECRs the 
matches are few: in part due to using a rather strict match 
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criterion of 2.5°. At heavy compositions a large fraction of 
UHECRs are matched to nearby extended radiogalaxies. 

For proton-UHECRs, of the seven models tested the BS- 
S and BS-A models maintain the high correlation between 
UHECRs and radiogalaxies in the Cen A region (Figs.[T]to[5]l. 
Accounting for deflections within the BS-S and BS-A mod- 
els lead to a higher concentration of UHECRs around Cen A 
(the UHECR detected toward Cen B is now also deflected up 
to Cen A) and the radiogalaxy WKK 4432. Applying deflec- 
tions from the models of SOS do not result in good matches 
between UHECRs and radiogalaxies in this region, though 
adding a dipole field to these models significantly improves 
the matches. In the case of the three matched events near 



NGC 7626, CGCG 413-019 and CGCG 514-050 (see Nagar 



& Matulich 2008 1, all of BS-S, BS-A, and AS-S are acceptable 



for UHECR-radiogalaxy matches; however, the SOS models - 
with or without a dipole component - do not provide consis- 
tently good matches. 

For many other UHECRs in Figures [T] to |5] we do not 
see obvious matches between the earth-arrival direction of 
UHECRs and nearby or intermediate (red and blue symbols) 
galaxies with radio jets. However, accounting for deflections 
by the GMF suggests several new potential matches between 
proton-UHECRs and extended radio jets. Some interesting po- 
tential matches between galaxies with distance < 200 Mpc and 
jet structures larger than ISO kpc are worth specifically not- 
ing: (a) we see possible matches between NGC 1275 and two 
AGASA-detected UHECRs; (b) CGCG 403-019, apart from 
a direct match to a R40-detected UHECR, is potentially the 
source of three additional AGASA-detected UHECRs; (c) the 
extragalactic Gamma-ray source with the highest maximum 
flux at HESS, Mi-k421, matches to an AGASA UHECR; (d) the 
BL Lacs H 2356-309 and lES 0347-121, both detected by 
HESS, can also be matched to AGASA events. Note that some 
of the AGASA events mentioned above do not appear in the 
figures as their 30% reduced energies are below 56 EeV. 

To explain the majority of UHECRs as originating in galax- 
ies with extended radio jets requires one of the following: (a) a 
stronger, or more ordered, intra-cluster and intergalactic mag- 
netic field resulting in extragalactic deflections greater than the 
current estimate (e.g. Harari et al. 2002 Dolag et al. 2004| l 



of a few degrees for UHECRs; (b) a larger mean free path 
or lower energy loss than current estimates of the GZK ef- 
fect, such that more distant radiogalaxies can contribute to the 
detected UHECR events; (c) a large fraction of heavy com- 
position UHECRs. The higher deflections of heavy UHECRs 
in Galactic and extragalactic magnetic fields result in differ- 
ences of several tens of degrees between the out-of-Galaxy and 
earth arrival directions. Recent results on UHECR composition 
from the Auger collaboration make this scenario likely (R\09a; 
PA09b). 

Given the small sample size of known UHECRs it is dif- 
ficult to evaluate how including UHECR deflections in the 
GMF eff'ects the correlation between UHECRs and nearby 



massive spiral galaxies (Ghisellini et al. 2008 1; this is best 



Gamma-Ray sources also revealed individually interesting co- 
incidences. However, their statistical significance is diflicult to 
evaluate and we do not discuss them here. 



4.6. Source Populations: Nearby Extended Radio 
Galaxies vs. Nearby Elliptical Galaxies 

Given the match between a subset of UHECRs and nearby ex- 
tended radio galaxies, it is relevant to test whether the potential 
association comes directly from the presence of extended ra- 
dio emission or from some underlying material associated with 
massive elliptical galaxies, for e.g. luminous or dark matter 
concentrations. As a preliminary test, we have used the RC3 
catalog ( [Corwin et al. 1994[ l to extract all nearby (D< 75 Mpc) 
elliptical galaxies. From this subset we identified all nearby 
RC3 ellipticals with extended radio jets and lobes, obtaining 
a list of 11 galaxies which includes all of our 15 nearby ex- 
tended radio galaxies except NGC1275, Cen B, WKK 4552, 
and CGCG 5 14-050. The limits in RC3 parameters for this sub- 
set of ellipticals are: distance D < 75 Mpc, morphological type 
-5 < T < -2, apparent B-band magnitude 7.69 < ma < 13.38, 
absolute B-band magnitude -22.6 < ms < -19.9, and size 
1.2 < D25 < 2.42, 0.02 < R25 < 0.2, 1.21 < Do < 2.45. We 
then selected a control sample of all RC3 ellipticals without 
extended radio jets and lobes which satisfied the above limits 
in RC3 measured parameters. This resulted in a control sample 
of 183 "radio-compact" ellipticals. Histograms comparing the 
distribution of distance, absolute magnitude, size, and angular 
separation from Cen A, of the two samples are shown in Fig.|6] 
The two samples are reasonably well matched with a small ten- 
dency for the radio-extended ellipticals to have slightly higher 
absolute magnitudes and to be slightly closer to Cen A on av- 
erage. 

A comparison of the statistics of matches of the two ellipti- 
cal samples to the outside-the-Galaxy arrival directions of PAO 
UHECRs is shown in Table|5] When normalized for the number 
of galaxies in each sample, the matches between UHECRs and 
extended radiogalaxies is significantly higher than in the case 
of radio-compact ellipticals. This diff'erence is most extreme in 
the BS-S and BS-A models. 



left for later analysis with a larger number of UHECRs. Our 
comparisons of UHECR outside-the-Galaxy arrival directions 
with the DRAGN sample, and with the several catalogs of 



4.7. Source Populations: Potential Galactic Plane 
Sources 

Very few UHECRs have been detected by PAO and AGASA 
in the region within 10° of the Galactic plane. Accounting for 
deflections in the GMF, however, results in several "Galactic 
plane crossing UHECRs" in most GMF models. In particular, 
the BS-S and AS-S models, and the models of SOS with or with- 
out a dipole component, serve to move the arrival directions of 
several UHECRs at -30 < b <0 towards the Galactic plane. 

We searched for potential matches between the UHECR 
trajectories within the Galaxy and objects from the catalogs of 
Galactic microquasars and Galactic magnetars. For the GMF 
models used, two or three UHECRs have trajectories which 
pass close (within 1 kpc) to a SGR or AXP: in the case of 
SGR 1900+14, a good match is obtained in all models ex- 
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cept the SOS models without a dipole component. In the case of 
SGR 1627-41, AS-S and AS-A model deflections lead to a con- 
nection to the mysterious close pair of events detected by PAO. 
Several UHECR trajectories also pass close to confirmed or 
candidate Galactic microquasars, and HESS-detected Galactic 
sources (e.g. Vela X and Vela Junior). Of course, fine tuning of 
the normalizations of the model permit closer matches; given 
the many free parameters (not least the uncertain distances to 
these Galactic sources, and the mass of the UHECR). We thus 
do not attempt a detailed study of the coincidences, but merely 
note that some UHECRs potentially cross the Galactic plane 
near candidate UHECR sources. 

5. Discussion & Concluding Remari<s 

The UHECR deflection simulations presented here, and the de- 
rived outside-the-Galaxy arrival directions of UHECRs, high- 
light several interesting points. For light (proton or Helium) 
UHECRs, the relatively small dispersions in the outside-the- 
Galaxy arrival directions for a specific UHECR and GMF 
model confirm that tantalizing insights into the source pop- 
ulation of UHECRs can be gained with a small sample of 
UHECRs. While heavy composition UHECRs suffer deflec- 
tions of several tens of degrees the concentration of outside- 
the-Galaxy arrival directions in the quadrant containing the 
SuperGalactic plane near Cen A also allows source population 
constraints. Conversely, picking one or a few UHECR source 
populations permits an evaluation of competing models of the 
ordered GMF. 

For light compositions, the trajectories of several UHECRs 
pass close to or through the Galactic plane within several GMF 
models. Interestingly, several pass relatively close to SGRs, 
AXPs, and/or microquasars in the Galaxy. |Ghisellini et al.| 
( 2008| l have argued that magnetars in nearby galaxies could be 
responsible for the PAO detected UHECRs. The integrated flux 
of HI emission from our galaxy is greater than the summed 
integrated HI fluxes of all massive spirals in the HIPASS sur- 
vey considered by lGhiselfini et al.| ( p008| l. If UHECRs origi- 
nate in magnetars in steady state - rather than in formation - 
then it would not be surprising to detect some magnetar-related 
UHECR events in the Galaxy. 



In Nagar & Matulich ( 2008 1 we claimed a correlation be- 



tween nearby galaxies with extended radio jets and the earth- 
arrival direction of a subset of UHECRs. By fortuitous chance 
these matched pairs lie along the lines of Galactic longitude 
which suff'er the minimum deflection by axisymmetric ordered 
Galactic magnetic fields. The correlation is strengthened with 
the results presented here. Accounting for deflections by BS- 
S and BS-A GMFs do not significantly change the previously 
claimed matches in the case of proton UHECRs. Rather, us- 
ing the BS-S model actually concentrates proton UHECRs 
closer to the nucleus of Cen A, and also results in several 
new UHECR-radiogalaxy matches. The models of SOS, with 
a dipole GMF component added, also maintain the previous 
matched pairs. Both the nearest radiogalaxy and the near- 
est radio-extended BL Lac are potentially sources of multiple 
UHECRs. It is possible but unlikely that all UHECRs origi- 
nate in Cen A. In the best case, i.e. choosing the most conve- 



nient global GMF model for any given match, about 30% of 
UHECRs outside the range -25° < b <0° could be matched to 
a galaxy with extended radio jets or an extragalactic Gamma- 
ray source in the case of proton only UHECRs. For heavy com- 
position or varied composition UHECRs we cannot rule out 
the possibility that all UHECRs originate in nearby galaxies 
with extended radio jets. We particularly noted several interest- 
ing and new matches between UHECRs and radiogalaxies for 
compositions near Oxygen. 

A remaining question is whether the correlation between 
radiogalaxies and UHECRs depends directly on the presence of 
extended radio jets and lobes or whether both trace an under- 
lying source population. We have briefly addressed this issue 
in a comparative test of nearby ellipticals with and without ex- 
tended radio structure: there is a strong indication that, within 
the BS-S and BS-A models of the GMF, the correlation with 
nearby extended radiogalaxies is directly related to the radio 
jets rather than to an underlying source population traced by 
massive elliptical galaxies. 

If radiogalaxies, or sources traced by radiogalaxies, are re- 
sponsible for a subset of UHECRs, a BS-S or BS-A model for 
the GMF is supported towards the Cen A region, i.e. the north- 
ern hemisphere of the fourth quadrant. The BS-S model (but 
not the BS-A) also allows a match between some UHECRs and 
Galactic plane sources. The models of SOS are consistent with 
both of the above only after we added a dipole field component 
(see Sect|2]for why this is dipole addition is not necessarily a 
self-consistent step). 

Finally, we remark that this work identifies the subset of 
UHECRs which can be matched to extended radiogalaxies us- 
ing a light composition (H or He) for UHECRs, and those 
which require a medium to heavy composition (e.g. Table [T]). 
The relative ratio of proton to heavy composition UHECRs 
then roughly agree with the composition mix suggested by 
PA09a and PA09b. It would be interesting to verify if the X^ax 
distributions of these two UHECR sub-samples confirm this 
composition difference. 
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Table 5 

Matches between PAO-UHECRs & Nearby Elliptical Samples 







RG Extended 


RC3 Extended 


RC3 Non-extended 


c 


Model 


Match 




Match 




Match 






BSS 


4 


7.4 


3 


9.4 


9 


1.2 




ASS 





0.3 





0.1 


12 


1.3 




BSA 


5 


7.1 


4 


8.8 


8 


1.2 


H 


ASA 


1 


0.1 


1 


0.1 


13 


1.4 




ASS+R 


1 


1.2 





0.1 


10 


1.3 




ASS+A 


2 


1.3 





0.1 


9 


1.2 




BSS(S) 


2 


1.7 





0.0 


9 


1.3 




BSS 


5 


5.6 


5 


7.3 


13 


1.8 




ASS 


2 


1.5 


1 


1.0 


9 


1.2 




BSA 


5 


7.5 


5 


9.5 


11 


1.6 


He 


ASA 


1 


0.5 


1 


0.7 


11 


1.5 




ASS+R 


5 


4.7 


2 


3.5 


6 


1.1 




ASS+A 


3 


2.3 





0.2 


5 


0.8 




BSS(S) 


1 


1.1 





0.3 


7 


0.9 




BSS 


10 


10.7 


10 


13.1 


18 


3.6 




ASS 


1 


5.1 


1 


6.1 


6 


2.4 




BSA 


7 


8.5 


7 


11.0 


15 


3.6 


O 


ASA 


1 


2.4 


1 


3.2 


5 


2.3 




ASS+R 


2 


5.9 


2 


5.5 


10 


2.6 




ASS+A 


3 


6.3 


3 


6.5 


9 


2.5 




BSS(S) 


2 


3.9 


2 


3.5 


10 


2.1 




BSS 


11 


13.3 


11 


17.0 


25 


4.8 




ASS 


2 


8.8 


1 


10.2 


9 


3.9 




BSA 


8 


10.7 


8 


13.3 


19 


4.9 


Al 


ASA 


1 


5.1 


1 


5.2 


10 


3.9 




ASS+R 


4 


8.7 


4 


9.2 


20 


3.7 




ASS+A 


1 


6.3 





5.6 


11 


3.7 









5.1 





3.5 


10 


3.0 




BSS 


15 


24.5 


14 


29.8 


19 


5.7 




ASS 


4 


16.2 


4 


19.6 


16 


4.8 




BSA 


12 


18.2 


12 


22.4 


20 


5.4 


Ca 


ASA 


2 


8.7 


1 


8.9 


10 


4.6 




ASS+R 


4 


9.7 


4 


11.3 


19 


4.6 




ASS+A 


3 


10.0 


3 


9.7 


16 


4.6 




BSS(S) 





8.1 





6.0 


13 


4.3 




BSS 


18 


31.2 


17 


37.8 


22 


5.9 




ASS 


7 


20.6 


7 


23.5 


19 


5.6 




BSA 


14 


24.9 


14 


29.2 


23 


5.9 


Fe 


ASA 


5 


15.9 


4 


16.1 


14 


5.5 




ASS+R 


6 


13.8 


6 


14.6 


23 


5.7 




ASS+A 


2 


12.3 


1 


11.2 


17 


5.3 




BSS(S) 


5 


13.1 


2 


10.1 


19 


5.3 



Note. — This table summarizes the match statistics between UHECRs de- 
tected by the Pierre Auger Observatory and three samples of ellipticals at D< 
75 Mpc: all 15 nearby extended (LLS > 180 kpc) radiogalaxies (columns 3 
and 4), all ellipticals from the RC3 catalog with extended (LLS > 180 kpc) 
radio structure (columns 5 and 6) and a matched sample of ellipticals from 
the RC3 catalog without extended radio structure (columns 7 and 8; see text). 
Columns 1 and 2 list the UHECR composition and GMF model. The two 
columns for each sample of ellipticals list the number of UHECRs (of the 
total of 27) which can be matched (see text) to some member of the sample, 
and the fraction of the total (1350) Monte Carlo points which fall within the 
match distance-limit of the sample's galaxies. 




Fig. 1. A comparison of the earth-arrival directions of all UHECRs with energies above 56 EeV detected by PAO (blue for 
UHECRs with energy > 75 EeV and green for UHECRs with energy between 56 and 75 EeV ) and AGASA (red open circles 
of radius 3.5°) with the estimated outside-the-Galaxy arrival directions of the same UHECRs (small colored dots) for our Monte 
Carlo simulations of the BS-S (red), AS-S (green), BS-A (blue), and AS-A (yellow) GMFs using a proton composition for the 
UHECRs. The colored circular symbols mark the positions of galaxies with radio jets at D < 75 Mpc (red), 75 Mpc < D < 
200 Mpc (blue), and 200 Mpc < D < 500 Mpc (black). In all redshift bins, filled circular symbols are used for galaxies with radio 
structures more extended than 180 kpc, and open circular symbols for galaxies with radio structures less extended than 180 kpc. 
Galactic SGRs (red triangles) and AXPs (blue triangles) are also plotted. The Supergalactic plane is marked by the dashed line. 
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Fig. 2. Same as Fig. [ijfor the three new models proposed by SOS: AS-S+RING (small red dots), AS-S+ARM (small green dots), 
and BS-S (small blue dots). Here, as in SOS, no ordered dipole field component was used. 
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Fig. 3. A comparison of the measured earth-arrival directions of all UHECRs with energies above 56 Ee V detected by PAO and 
AGASA (blue, green, and red open circles with radius 3.5° as in Fig{T]i with the estimated outside-the-Galaxy arrival directions 
of the same UHECR (crosses) within the BS-S model for UHECR compositions of protons (red). He (green). Oxygen (blue). 
Aluminum (yellow), and Calcium and Iron (black). The first three compositions of the same UHECR are connected with blue 
lines while the heavier compositions are not connected to avoid overcrowding. Other symbols and names in the figure are the 
same as in Fig.[T] 
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Fig. 4. Same as Fig. [l]but for a UHECR composition of Oxygen. 



Nagar & Matulich: UHECRs: Galactic Magnetic Field Deflections 



21 





20 40 60 80 
Dist. (Mpc) 





70 








60 


' 1 






50 


' 1 








1 , 


w 


O 






o 




40 












? 


O 






o 




30 






d 






d 


Z 






z 




20 


r ~ 






10 


' ""i - 














-23 -22 -21 
Mb (mag) 



-20 



40 



20 



10 



12 3 
D(CenA) (rad.) 



Fig. 6. Histograms of Morphological Type, Distance, Absolute B magnitude and distance from Cen A for the two samples of D 
< 75 Mpc elliptical galaxies drawn from the RC3 catalog (see text): the 1 1 nearby RC3 ellipticals with extended radio structures 
(solid line histograms) and the matched control sample of 183 nearby RC3 ellipticals without extended radio structures. 



